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Bv Percy J. Bobbitt 


SUMMARY 

The lifting-surface sidewash behind rolling triangular wings 
has been derived for a range of supersonic Mach numbers for 
which the wing leading edges remain swept behind the Mach 
cone emanating from the wing apex. Variations of the side- 
wash with longitudinal distance in the vertical plane of sym- 
metry are presented in graphical form. 

An approximate expression for the sidewash has been devel- 
oped by means of an approach using a horseshoe-vortex approxi- 
mate-lifting-line theory. By use of this approximate expres- 
sion , sidewash may be computed for wings of arbitrary plan 
form and span loading. A comparison of the sidewash com- 
puted by lifting-surface and lifting-line expressions for the 
triangular wing shoived good agreement except in the vicinity of 
the trailing edge when the leading edge approached the sonic 
condition. 

An illustrative calculation has been made of the force induced 
by the wing sidewash on a vertical tail located in various longi- 
tudinal positions. 

INTRODUCTION 

In order to make reliable estimates of the total forces and 
moments acting on an aircraft, accurate evaluations arc 
required of the loadings on the individual isolated compo- 
nents and of the interference effects between components. 
Although considerable effort lias been expended in recent 
years to supply much of this needed information for the 
supersonic speed range, many important problems remain. 
Among these is the induced effect of the wing flow field or, 
more precisely, the wing sidewash on the vertical tail. The 
only specific numerical results of this nature obtained to 
date have been for the angle-of-attack motion. In the 
vertical plane of symmetry for this case, however, the side- 
wash is zero and tail surfaces located in this plane are un- 
affected. This is not the situation for the rolling, yawing, 
and sideslipping motions where the sidewash in the vertical 
plane of symmetry is finite and the load induced on the 
vertical tail can be appreciable. Evaluation of the sidewash 
for these motions would, therefore, be important in the 
prediction of the lateral stability of supersonic aircraft. 

The present report presents the derivation of the sidewash 
behind steady rolling, triangular wings with subsonic leading 
edges. Both lifting-surface and lifting-line methods, pre- 
viously applied primarily to determine downwash, are utilized 


and comparisons are made of the sidewash computed by the 
two methods in order to give an indication of the worth of the 
more easily obtainable lifting-line results. The lifting- 
surface sidewash is determined by using the doublet-distri- 
bution method of reference 1, and the lifting-line values 
are obtained by use of the lifting-line approach given in 
reference 2. 

An illustrative calculation using the derived sidewash is 
made of the force induced on a half-delta tail operating 
behind a rolling triangular wing, and this force is contrasted 
to the force that would act on the tail if it were rolling in 
the undisturbed stream. 

The material presented in this report was submitted to 
the University of Virginia as a thesis in partial fulfillment of 
the requirements for the degree of Master of Science in 
Aeronautical Engineering. 


SYMBOLS 


The positive directions of forces, moments, and velocities 
are shown in figure 1. 

Cartesian coordinates of field point 
Cartesian coordinates of doublet or line-vortex 
position 

perturbation velocities along x-, y -, and 2 -axis, 
respectively 
wing aspect ratio, b 2 /S 
wing span 

Yawing moment 


z, V , 2 

2/i> 

U , V , w 


A 

b 

C n 


yawing-moment coefficient, 


gSb 


Gn r— 

C Y 

Cy = 

c 

d 

h 

hi, h 2 


dCn 
2 V 


dCy 

d P b 

2V 


p — »0 


side-force coefficient, 


Side force 


p — K) 


wing root chord 

distance from wing trailing edge to a point 
downstream 

displacement of vortex sheet below wing trailing 


edge 


limits of ^-integration 


1 Supersedes NACA Technical Note 3009 by Percy .T. Bobbitt, 1950. 
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Figure 1 . —Triangular wing oriented with respect to body system of 
axes used in analysis and associated symbol data. 

i variable index used in summations and as 

subscript 

k constant 


v / (x-*i) 2 -ff 2 z 2 

00^1 

V(x— c) 2 — P 2 Z 2 


OoC 

ki,o'=\ 1 — &i,o 2 


y, v 7 sidewosh 


Yawing moment 


k 2 -~ 


OoX i 


■ylix-xtf-pZ 2 

OoC 


k 2 , o — 


^(x-c) 2 -p 2 z 2 


k 2t o' — VI — & 2 , ()" 

M free-stream Mach number, V7 Velocity of sound 

in free stream 

m slope of lifting line 

P static pressure 

AP=P i — P u 

p angular velocity of roll, radians/sec 

q free-stream dynamic pressure, ^pl 72 

S wing area 

A u s =n u —Ui 

V free-stream velocity 

v P side wash induced by doublets distributed over 

plan form 


Vw 


X=x—X\ 

Xi=x—Xi 

x-X± 

‘~b/2 

X 

x,=- 

r - Xl 

Y=y—yi 

Yi=y—yi 

Y -Zi 

6/2 


sidewasli induced by doublets distributed over 
wake 


2/o= 

Zo= 


y 

6/2 

z 

6/2 


a angle of attack, radians 

fl=^M 2 - 1 

T circulation at any spanwise station 

e angle of downwasli measured in a; 2 -plane, be- 

tween trailing vortex sheet and axis parallel 
to free-stream direction, radians 


=@ f° r t-riangular wings^ 


00 

c \ 

X angle through which vortex sheet rotates in 

moving from wing trailing edge to a point 
which is cl distance downstream 
p density of free-stream air 

<#> perturbation velocity potential 

A </>$ 4>u 

= 0o_c =s I 

fiz Zo 


E 


complete elliptic integral of second kind with 
modulus k , 


, rv-fcv 

"Jo vr= 


ds 


E(t,k) 
F(t,k ) 


complete elliptic integral of first kind with modu- 

llls *J. vi-ivyi-s 2 

incomplete elliptic integral of second kind with 
argument t and modulus k , [ . —f S n - , -■ 

Jo ^li-k 2 s 2 ^l-s 2 


incomplete elliptic integral of first kind with 
argument t and modulus k, 


, C l -y 1 —k 2 s 2 ds 
Jo -y/ 1 — 


, x—sl6 0 2 x' 2 +P' 2 z 2 (l—do 2 ) 

i-e 0 2 


f 0= 


J'n — 0t>\ rp' d~ 2 q* ( 1 ft)~ ) 


“^ = 9 _ l _ (^ 2,0 — E 2 ,o)F^COt 1 — ) k 2 , o') 

# 2 , 0 # (cot- 1 k 2 , o'/ 

2 — 0 n 2 A- 2 

G%): 


e 0 2 


i-0„ 2 E, ^~r=e? F ' {9o) 
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H=-g- 

P 2 G(e o) 

L.E. leading edge 

T.E. trailing edge 

| sign denoting finite part of integral 


Subscripts: 


D 

E 

l 

P 

s 

u 

W 


conditions in region D (fig. 2) 
conditions in region E (fig. 2) 
pertaining to lower side of surface 
plan form 

conditions on surface of discontinuity (at z x = 0) 

pertaining to upper side of surface 

wake 


The subscripts 1,2, 1,0, and 2,0 on the elliptic functions E 
and K indicate the modulus of the elliptic function ; that is, 


*,) 



ANALYSIS 

GENERAL REMARKS 

The problem to be considered herein is that of determining 
the perturbation sidewash velocity behind a rolling triangular 
wing for a range of supersonic Mach numbers for which the 
leading edges of the wing are subsonic. The analysis is 
based on an application of linearized supersonic-flow theory 
and, hence, the results obtained will be valid within the 
limitations of linear theory. 

In the analysis which follows, several assumptions are 
made concerning the trailing vortex sheet. These assump- 
tions are that the vortex sheet must remain flat- behind the 
wing and that the rotation of the vortex sheet is small enough 
to be neglected. In addition, the nonrestrictive stipulation 
is made that the rolling wings be at zero angle of attack. 
Further discussion of these points will be found in the section 
entitled “Results and Discussion.” 

In flight, a steady rolling motion will usually be maintained 
by differentially deflected ailerons that create a sidewash 
opposed to the wing sidewash. Calculation of aileron side- 
wash, which may be of the same order of magnitude as the 
wing sidewash, will not be considered in the present report. 

BOUNDARY CONDITIONS 

The boundary conditions for the proposed problem may 
be prescribed on the z = 0 plane and are similar to those given 
for the angle-of-attack motion in reference 1 . 

The downwasli boundary condition on the rolling wing is 


Pressures on the wing and pressure differences across the 
wing surface are finite and, for a great variety of plan 
forms, have already been obtained. (See, for example, 
refs. 3 and 4.) Off the wing and in the plane of the wing, 
the pressure, and hence the pressure difference, must be zero. 

In the z = 0 plane, the local pressure difference is directly 
proportional to the streamwise component of the perturba- 
tion velocity and is given simply as 

AP{x u y x ) = 2 Au s (x u y x ) /jx 

<i v 

By consideration of the relationship between the perturba- 
tion velocity potential and the streamwise velocity com- 
ponent; that is, 


u,(x u yi) ( 2 ) 

an expression giving the jump in velocity potential across the 
# 2 /-plane in terms of the local pressure difference may be 
written as 


A0,(si,2/i) 


aP(x, 
2 J l.e. #o 


Vi) 


dx 


(3) 


Since, from equations (1) and (2), 

AP(x u y l )_2 DAhixiffli) m 

q ~V K } 

and since &(#i,^i) is an odd function in z , it is clear that, 
beyond the trailing edge, A0 S must be independent of x x 
to satisfy the zero-pressure condition in the wake. The 
integration indicated in equation (3) should, therefore, be 
made from the wing leading edge to the trailing edge to 
obtain A <t> s in the wake. 

SOLUTION TO BOUNDARY-VALUE PROBLEM 

The linearized partial-differential equation which the per- 
turbation velocity potential must satisfy in supersonic flow is 


2 b 2 4> d 2 0 d 2 <t> 

P dx 2 by 2 bz 2 ~ 


(5) 


For the problem being considered herein, the solution to 
equation (5) may be written as 


<Kx,y,z)= 


zj?_ CT A<t> s (x\,y x )dx\ dy x 

2tt JJ [(x—XxY—P^y-ytf-pz 2 ]* 12 

Plan form 


w= (py)p^o 

In order to analyze the quasi-steady rolling problem by use 
of steadv-flow theory, the rolling wing is considered fixed in 
approximately the 3=0 plane but twisted linearly in the 
spanwise direction. Only small linear twists are allowable, 
however, in order not to violate the assumptions of small- 
perturbation linearized theory; hence, the rate of roll is 
necessarily small (approaching zero). 


zj- CC A4>Ayi)dx, d i ji 

2w JJ [{x-xX-pty-yX-reY* 

Wake 


( 6 ) 


This expression represents the potential in space due to a 
distribution of doublets in the xt/-plane with strengths that 
are governed by the potential jump across the 3=0 plane. 

The symbol | indicates that the finite parts of the 

infinite integrals are to be taken when they appear. 
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APPLICATION TO TRIANGULAR WINGS 


The loading over a rolling delta wing with subsonic lead- 
ing edges has been found in reference 4 to be 


AP, ^ x 4 ^100-/3 ?/! 

2 J Vp 2 G(d 0 ) V*o V-0 V 

where 

and 

F'{do)=F(^,^e^ 


(7) 


From equations (3) and (7) the potential jump across the 
wing surface is 

A <t> s (x l ,y l )=H(3iji ■\Q* 2 x l 2 —& 2 y* (8) 


( d<f>w\ 

Vw WyyU 


zHp 2 5 

f f 

13 Vi -J6 2 c 2 — fi 2 y 2 dx x dy x 

2ir dy 

JJ 

Wake 

[(x—x x ) 2 —fi 2 (y—y x ) 2 —fi 2 z 2 } 3/2 


Subsequently it will be convenient for computational 
purposes to derive expressions for the sidewash which have 
been nondimensionalized by pb/ 2 so that 


v v P ■ Vw 
pb/ 2 pb/2'pb/2 

When ! is written in a slightly different, though equiva- 
pb/2 

lent, form 

v/V 

pb/2V 


and in the wake 


where 


A<fc( 2 /i)=# 02 /i a d 0 2 c 2 —P 2 yi 2 


H= 2p 
P 2 G(d 0 ) 


(9) 

( 10 ) 


The velocity potential in space may now be written as the 
sum of the two expressions (see eq. (6)) 


_zHp 2 

rr Pyi Ydo 2 x x 2 ~l3 2 y x 2 dx x dy x 

2t 

J J [(x-x x ) 2 -p 2 (y-y x ) 2 -(3 2 zT 2 

Plan form 


</> P = — 


and 


( 11 ) 


<t>w= 


zHp 2 r r Py\^0o 2 c 2 —P 2 yi 2 dx x dyi 
2tt J J [(x— x x ) 2 — /3 2 (t/ — yd 2 — p 2 z 2 ] 312 

Wake 


( 12 ) 


As previously stated, the primary purpose of this report is 
the determination of the velocity perturbated behind the 
wing parallel to the y-axis (or the sidewash). This flow 
velocity may be obtained by taking the partial derivative of 
the velocity potential with respect to y, or 


With 


dy 




it can be recognized that the nondimensional sidewash 
parameter may be defined as the induced angle of sidewash 
per unit wing-tip helix angle pb/2V. 

The rest of this section is devoted to the evaluation of 
equations (13) and (14) at points in the following two 
regions of the x^-plane (see fig. 2): 

(1) The region lying between the Mach lines emanating 
from the wing trailing edge and the line of intersection of 
the two cones from the trailing-edge tips. 

(2) The region which extends from the line of intersec- 
tion of the two cones from the trailing-edge tips downstream 
to infinity. 

These two regions are denoted, as in reference 1, by E and 
D, respectively. The contributions of the doublets dis- 
tributed over the plan form and the wake to the sidewash 
in regions E and D are considered separately. 

Sidewash due to doublets distributed on plan form in 
region D. — In region 1) the sidewash contributed by the 
doublets distributed on the plan form is 

d<t) PtD r ( zH& 2 d I' f 0 Py\ V flo W—P 2 yi 2 dyi dx l \ 

dy 2tt Jo J e 0 x l [(x—x l ) 2 —p 2 (y—y 1 ) 2 —p 2 #] z/2 J 

(15) 

In order to facilitate the integrations involved in deter- 
mining v PfD , it is convenient to carry out the differentiation 
and limiting processes first. This procedure gives 


the sidewash in the xz-plane from equations (11) and (12) 
will be given by the sum of 


v P = 


/ d 4) A 
V d y A— >0 


=— lim 

2/ — *0 


zHp 2 d^ rr Pyi^QoW—Fyfdx! dyi 
2?r dy J J [(x—x 1 ) 2 -p 2 (y-yi) 2 —p 2 z 2 ] z/2 

Plan form 


(13) 


3 zHfi 2 

Vp, D — * 

7 r 





e P 3 Vi 2 yd 0 2 x x 2 —P 2 y x 2 dy x 

[{x-xtf-pyS-pzT 2 


(16) 


When the following substitutions are made, 

, Fvl 

8 ~e 0 W 

, t e 0 % 2 

2 (x-Xif-pz 2 
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where 




OoW 


0 o 2 X it o 2 


2 (x—X^—^Z 2 (Xo—X lt o) 2 —do^o 2 


Sidewash due to doublets distributed on plan form in 
region E. The portion of the wing area over which the 
integration in equation (13) is to be performed is different 
for each position of the field point in region E. This fact is 
evidenced by the appearance of the field-point coordinates 
x , y, and z in the limits of integration. It is expedient in 
determining v PfE to follow the same procedure used in deter- 
mining VP' D (differentiating before integrating). This is 
allowable since it can be shown that the expression for the 
potential <t> PtE can be differentiated with respect to y without 
regard to the variable limits when the evaluation of the 
derivative is made at y= 0. 

Differentiating <t> PiE with respect to y and then setting y 
equal to zero yields 


v P , s'- 


*** (w r 

7T Jo Jo 


owl 2 M^r 

■x y^-p- 


-yNy\ 


r(x—xi) 2 —p 2 z 2 2 t /2 

L y ' J 


Line of intersection of 



H c + 


Figure 2. — Regions behind a triangular wing. 


v P , D becomes 


3 zHfi 2 r c k 2 5 dxi C 1 s 2 VI —s 2 ds 

7T Jo BqXi Jo (1 —k 2 2 s 2 ) 5/2 


(17) 


which, by the use of partial fractions, may be written in the 
more amenable form 


3 zHp 2 r k^dxi r r ds 
Vp - D ~ T Jo e oXl L Jo VT^VT 31 ^ 2 

n ds _ _ 

{2 ~ h2) J a ^ 1— <i*(l —k 2 2 s 2 ) 3,2 ~ 

(i -fe 2 )j o ~ 1 _ s 2( 1 _^ 2 2 s 2)5/ 2 ] (18) 

The integrals in equation (18) may be reduced to standard 
elliptic forms by use of the Jacobian transformation, s=sn u 
(refs. 5 and 6), and readily integrated to give 




dx i 


(19) 


By replacing H by its equivalent and nondimensionalizing 
Xi, x, and z, equation (19) becomes 


r p, p 2gp f* k’2 / 

pb 2~~G(0 o)tJo 0oXi,o \ 


-2 1 < 2 + 


2 ~^ 2 E J 

1 -fc/v 


dxi.o 


(20) 


3 zHp 

7 T 




tJ(X-Xi)*-PZ* y2 


P 


, / OoW 

v F 


-yddyi 


nx-xi) 2 —p 2 z 2 2 t /2 

L p y> J 


( 21 ) 


where 


f= 


-tJ0 o 2 x 2 +I3 2 z 2 ( 1 - 0 „ 2 ) 
l-0o 2 


The y r integration in the first double integral of equation 
(21) is identical to the ^-integration in equation (16); hence, 
only the second term of equation (21) remains to be con- 
sidered. 

In the integration 


I 


let 


where 


y>^W_ y 2 dy 




^ k 2 t 2 =y 2 


k i 2 = 


(.x—XiY—i 8 2 z 2 


0o 2 *i 2 


This substitution results in the expression 


0 I f 1 t 2 -\/l — k 2 t 2 dt 
0oX\k\ 2 Jo (l-< 2 ) 5/2 


( 22 ) 


which, except for the use of the finite-part concept, may be 
integrated in a manner similar to that used for equation (17). 
Performing the integration in expression (22) allows v P , E 
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(eq. (21)) to be written as 

zHi 3 2 P , / jr ,2-k-f „ \ 

B= ^TJo k 2 (“ 2 ^+ 1 ^J 4 

2/ZpPrf/, 1 ,1-2*,* \ 

- J, 1 + l -*, 2 W 


(23) 


Equation (23) for v P% E does not lend itself readily to numerical 
calculations because, at the point x x =f, k x and k 2 become equal 
to 1 and give rise to first-order infinities in the integrands. 
Appendix A shows how this difficulty is eliminated bv a 
parts integration of the terms containing the singularities. 
After the singularities have been treated and the x x -,x- } and 
2-lengths nondimensionalized, equation (23) becomes 


v P , 


pb/ 2 


_ 2z 0 f Cfodx^ k( _ 2K , F) , ndx t , o -K l +E 1 _ p 0 

'*<?«?„) \ Jo ^.o 2( 2+ 2 )+ J/o^..o h 2 Jo 


— r -q-vk.+*i 

L^A!±2o!(1z^J 

do V*o 2 +2o 2 (1— 0?) 


[' 


t tnb -r-( i -«.‘K.±&l 

MTifa) - *i.oTI p I ^-vV+go^l-tfoOj f t - 


k 2 ^O-Tl.O 

t ^x£-0£z£ _ coth _, 




o~ dp^Zp' 
2<9oA'V+2 0 2 (1— 0o 2 ) 




#o ^ Xq 2 -\- Zq( \ — 0 O 2 ) 
I Ei o 


0 O V ^o 2 + ^o 2 ( 1 — 0o 2 ) V X 0 2 + Zo 2 ( 1 — 0o : 


[ ka ~ - 

2 ) H~ x 0 

dp^Xp 2 -\-Zp 2 (\ d 0 2 ) j 


tanh-‘ 


(24) 


where 


X __ Xo — do VV + 2 0 2 ( 1 — do 2 ) 
h i -d 0 2 


Numerical calculations may be made by using equation (24) for all values of 0 O except 0 O = 1 where / 0 becomes indeterminate 
and the arc hyperbolic functions become infinite. The indeterminacy when evaluated yields 


/ X \ xo 2 -doW 
Uo) e o=i 2x 0 

Thus, an integration by parts of the singular terms of equation (23) similar to that made to obtain equation (24), using now 
(Jo)e 0 =i and expressions for k x and k 2 in which 0 O has been set equal to 1, yields 


Vp, 

pb/2 




Ox" J f (2*0^1. o — #o 2 +Zo 2 ) j 


(E x —K x )(xo—x i,o)' 


k 2 x Xt o 


] 


dx, o- E ^ kl ~^° — — log, ( 2 r 0 -x 0 2 + 2 0 2 )+ 

2xo 


7i-y.ro 2 — 2o 2 

4x 0 


log, (x 0 : 


r-2o 2 )} 


(25) 


The integrations in equation (25) may be handled by numerical methods. 

Sidewash due to doublets distributed in wake in region D. In region D, equation (14) takes the form 


v W' D — lim ~ 
1/ — *o 


zH& 2 d 


2 x dy 


'**-0V(tf-I'l) 2 +Z 2 


dx i 


[(.x-XtY-piy-ytY-Fz T 2 


f PyiyOoV—pyfdyi f 

J _0OC J c 

~T 

Carrying out the differentiation and then the first integration for ;/=() gives 

(3-c) 2 ffWVfloV-|8W 


(26) 


V\V,D — “ 


zH(x-c) r e 


7T/3 2 





(y 2 +z 2 ) [(x— c) 2 — /3v/i 2 — p 2 2 ] 1/2 (yi 2 +z 2 y [(x— c) 2 — /3 2 t/! 2 — p 2 z 2 ] 3 ' 2 


~d(Pyi) 


(27) 


By making the variable substitution 


and defining 


Vi 


e 0 2 c 2 


e 0 V 


(x— c) 2 — 13 2 2 2 
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equation (27) reduces to 
c)9 0 3 c 3 ^2. q f 1 


V\Y, D~ 


w P 2 Z* 
H(x—c) 3 0ack 2 'ii 


' f 1 s 2 yi 

Jo / 1+ *oV 

V + 2 2 /3 2 


-s 2 ds 


^ \ 1 — & 2 . OS 2 
-cj%>chX f 1 s 2 V I — s 2 ds 


(28) 


The expression for iv,z> given by equation (28) may be 
integrated (see appendix B) to yield a closed-form solution 

for v WiD . The nondimensional equation for is 


Pb / 2 


V\V,D_ 


pb/2 rG(d o) 


^ ^+(K 2 .o-E 2 . o) F (cot- 1 1, fc 2 .o' )- 


K^oEfcot- 1 *, 0 'Y1 °Xf °<f ° 1)z ° \ 

V *o ’ /JV^+1 ^o(2o 2 +1) / 


/ 


(29) 


Sidewash due to doublets distributed in wake of region 

E. In region E the derivation of the sidewash is similar to 
that of region D and yields 


Vw, e 2 (#q 1 ) 

pb/2 7r0o6?(0o) 


{[i- 


K\, oflo^oV £q 2 H~ 1 
x 0 — 1 


(Ki. o-E li0 )F (cot- 1 *!.#')- 
tf,.oE ( cot- 1 - (22 ° 2+1)e ° 2K 1 .oZo\ 

V So /JVz 0 2 +l (* 0 - 1 ) / 


(30) 


where 


ho 2 


(x—c) 2 —fi 2 z 2 

0o 2 c 2 


Sidewash at a>axis. — In the xy- plane (plane of the wake), 
only the doublets distributed in the wake contribute to the 

sidewash. Equation (20) for and equations (24) and (25) 

for approach zero as z 0 approaches zero, whereas equa- 
tions (29) and (30) give 


Vw, d Vw, E 


pb/2 pb/2 6(6 0 ) 


(31) 


This result is identical with that which would be obtained 
by use of the formula obtained in reference 7 by considering 
the properties of vortex sheets. This formula, for y=y l = 0, 


pb/2 


(*, 0 , 0 +) 


— ■ 


pb \dyj v .. o 
21 2~ 


where 


r — ( a</>, s )te 


(32) 


(33) 


Sidewash at x= oo (Trefftz plane). — As x approaches 
infinity, the contribution of the doublets distributed on the 
plan form to the sidewash goes to zero and the total side- 
wash is given by 

418784—57 2 


(34> 

Equation (34) could also have been determined more di- 
rectly by using the formula 


( v \ ;• g |~d f 6/2 A^jxr^y^dyn 

\P*/2A-.» 2wpbj2 |_dy J_>,/2 (y—y\Y+z 2 J 


(35) 


which may be obtained from equation (26) by performing 
the first integration and then setting x equal to infinity. 


LIFTING-LINE SIDEWASH 

The lifting-surface method by which the sidewash behind 
a rolling delta wing was derived in the previous section is 
applicable to wings of arbitrary plan form; however, the 
integrations which would be required before the potential 
or one of the perturbation velocities could be obtained in a 
calculable form are extremely difficult to evaluate. It is of 
importance, therefore, to develop some approximate ex- 
pressions which may be easily evaluated either analytically 
or numerically. References 2 and 8 indicate that a lifting 
line and a lifting line approximated by supersonic horseshoe 
vortices can be used as good approximations to lifting- 
surface solutions for most downwash problems. It would 
seem that a comparison of the sidewash behind the rolling 
delta wing calculated by the lifting-surface method with that 
calculated by an approximate method might give some 
indication of the usefulness of the “approximate” approach 
for sidewash problems. 

References 2 and 8 together represent a fairly thorough 
study of the lifting-line and approximate lifting-line methods, 
especially with regard to downwash calculations, and show 
that the swept lifting lines will probably give the best 
results for swept and triangular wings. The potential due 
to a yawed lifting line may be obtained from the errata of 
reference 2 as 


nih) z^X 2 -p 2 (Y 2 +z 2 ) 

2t Y X — - — — 

m m 


J_ p2 dV(yd 

2t Jh , <&/, 


tan 1 


zYX 2 -p 2 (Y 2 +z 2 ) 



XI 

m 


dy i 


(36) 


where the equation of the lifting line is 

_Vi±k 


Xi'- 


an 


and the circulation r is defined by equation (33) as the 
potential jump across the surface evaluated at the trailing 

edge. (See fig. 3.) When is zero, equation (36) 


dy i 


becomes 


YX-tXl 

m m 


2i r 


( 37 ) 
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which might be considered as the potential in space (at a point x, y, z ) of a linite yawed vortex of constant strength. 
A number of finite yawed vortices distributed along a line can be used to approximate the potential in space of a lifting 
line with any prescribed lift distribution. 

From equation (37), the sidewash due to a yawed vortex is readily obtained by taking the derivative with respect to y. 
The following result is obtained: 


dy 2 tt 


zY_ 

m 


(2 X' 2 -p 2 Y 2 -p 2 z 2 )-zX(X 2 -l3 2 z 2 ) 


JX 2 -pHY 2 +z 2 )[(YX-~- Y ^y+z 2 (X 2 -p 2 Y*-l3 2 z 2 )J 


(38) 


When m approaches infinity, equation (38) becomes the sidewash for a rectangular horseshoe vortex and agrees with the 
equations given in references 9 and 10. 

Since the loading on a rolling wing is an Asymmetrical, the induced sidewash from each panel is in the same direction 
and equal in the y=0 plane. For this reason it is necessary to calculate the sidewash only from one panel and double it. 

Equation (38) can be utilized to formulate an approximate expression for the sidewash due to a series of constant-strength 
yawed horseshoe vortices spaced along a line so as to represent as closely as possible the span load distribution due to rolling. 
This expression is 


^ r(y <+1 )-r(y^ t ) 

1 U 4ir 


zY 

m 


H2X ( 2 -0 2 Y i 2 -P 2 z 2 )-zX ( (X i 2 -0 2 z 2 ) 


^-^-^z^Y.X.-^-^y+z^X, 2 -^ 2 -^ 2 )^ 


(39) 


where Y t =y—y it X t =x—x ( , and subscript i takes on all integral values from 0 to n. Equation (39) in non- 

dimensional form is 


pb/2 




,=u imjxs-p 2 Yi, 0 2 -f3W (X^-^Y^-ph/)] 


(40) 


In the application of equation (40) to the calculation of sidewash, some care should be exercised that the forecone from 
the field point under consideration does not intersect the lifting line at a point close to the corner of a yawed horseshoe vortex. 
When the forecone intersects the lifting line near a corner located within the forecone, the expression under the radical in 
the denominator of equation (40) becomes small and the sidewash becomes large. (See sketch 1.) A zero value for the 
square root and an infinite value for the sidewash result when the forecone intersects the corner. The abruptness of the 
infinity varies with the distance of the field point from the corner. Note that when x= o° , the infinity no longer exists. 
The preferable field-point locations have forecones intersecting the lifting line as shown in sketch 2. The closeness of the 
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( 0 , 0 ) 


y>y i 



Figure 3. -Finite yawed vortex used to approximate a swept lifting line. 

forecone to the corner when the corner is outside the Mach 
cone is obviously of no consequence because only line 
vortices within the forecone from the field point contribute 
to the sidewash at the point. 


RESULTS AND DISCUSSION 


EXACT SIDEWASH 


The exact linearized nondimensional sidewash has been cal- 
culated for values of 0 o of 1.00, 0.75, 0.50, 0.40, and 0.30; values 
of^o from 1.2 to 2.4; and values of z 0 from 0 to 0.6 except 
where these values are ahead of region E. Variations of the 


sidewash parameter 


v/V 


with # 0 for 7 values of z 0 from 0 to 


pb/2V 

0.6 and for the d Q ’s given are presented in figure 4. Cross 
plots of figure 4 which show the variation of the sidewash 
parameter with z 0 for 7 values of x 0 from 1.2 to 2.4 are given 
as figure 5. 

In order to depict the effect of Mach number and leading- 
edge sweep, variations of the sidewash parameter with z 0 
for values of 0 O of 1.00, 0.75, 0.50, and 0.30 have been plotted 
for three longitudinal locations: # 0 — 1-6, a: 0 =2.0, and x 0 = c° 
(fig. 6). (It should be noted at this point that an increase 
in 0 O may be interpreted as either an increase in Mach num- 
ber for a fixed leading-edge slope or an increase in the wing 
semiapex angle for a specific Mach number.) The major 
difference to be noted in the effects of changing 6 0 is that, 
when the longitudinal station is ahead of the line of inter- 
section of the Mach cones from the trailing-edge tips, an 
increase in 0 O causes an increase in the sidewash at the higher 
values of z 0 , whereas the sidewash at a station remaining 
behind the intersection line during an increase in 0 O experi- 
ences a decrease in sidewash at all values of z 0 which are un- 
affected by the localized infinity at the intersection line. 
By way of illustration it can be seen at station #0=1.6 that 



(b) 0„=O.75. 

Figure 4. — Longitudinal variation of nondimensional sidewash param- 
eter in XoZo-p\ane behind triangular wing for a number of values 

pb/2V 

of zq and 0„. 

when 0 O is increased from 0.75 to 1.00 the sidewash increases 
for values of z 0 greater than 0.2; for values of 0 O of 0.3 and 
0.5, when # 0 =l-6 is behind the intersection line, the effect 
of increasing 0 O is to decrease the sidewash at all values of 
z 0 except at z 0 =0.6. This point is affected by the infinity at 
the intersection line. 

APPROXIMATE SIDEWASH 

From the nature of the analytical and numerical integra- 
tions required to obtain the exact sidewash for triangular 
wings, it is apparent that for wings with more complex 
potential-jump expressions the derivation of exact sidewash 
would be a difficult task. Herein lies the merit of the ap- 
proximate lifting-line method (eq. (40)) which is not en- 
cumbered by the complexity of the wing-loading expression. 
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(c) 0o=O.5O. 

(d) 0o=O.4O. 
Figure 4. — Continued. 


Nondimensional sidewash parameter, 
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(a) 0o=l.OO. 

Figure 5. — Variation with z 0 of nondimensional sidewash parameter 
in £ 0 Zo-plane at a number of longitudinal stations. 



0 


.1 .2 



.5 .6 


(b) 0o=O.75. 
Figure 5. — Continued. 


Nondimensionol sidewash 
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(c) 0o=O.5O. 
Figure 5. — Continued. 


(d) 0o=O.4O. 
Figure 5. — Continued. 


Nondimensional sidewash parameter, 
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(e) 0u=O.3O. 
Figure 5. — Concluded. 



(a) x Q = 1.6. 

Figure 6. — Change in vertical variation of sidewash with 0<> at x c = 1.6, 
2.0, and infinity. 



(b) *o = 2.0. 

Figure 6. — Continued. 
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(c) x*a=<». 

Figure 6. — Concluded. 

The approximate method is, however, hindered to the extent 
that an area distribution of loading is assumed concentrated 
on one or several lines. The penalty that this assumption 
imposes on the quality of the results cannot be ascertained 
in every case. For the triangular wing treated herein, com- 
parisons may be made between the results from the lifting- 
line and the lifting-surface methods, and perhaps some indi- 
cation may be obtained as to the regions wherein the approx- 
imate method may or may not give reliable values. 

By use of equation (40), approximate lifting-line calcula- 
tions of the sidewasli have been made for 0 o =l-OO and 0.40, 
values of x 0 from 1.2 to 2.4, and values of z 0 from 0.1 to 0.6. 
Sidewasli values for z 0 =0 were obtained from equation (32). 
A comparison of the sidewasli calculated by the lifting-line 
and lifting-surface methods has been made in figure 7, and 
the agreement is shown to be good everywhere except at the 
high values of z 0 close behind the trailing edge for 00=1 TO. 
The agreement in this region is considerably better at 
0 o =O.4O (fig. 7(b)) and indicates that as 0 O is decreased 
from 1.00 the approximate calculations will become more 
reliable at locations close behind the trailing edge. 

Seventeen yawed' horseshoe vortices were used to approxi- 
mate the lifting line, with the concentration of vortices 
greater near the tip because of the rapid change in the span 
loading in this region. The lifting line used in the ap- 
proximate calculations consisted of a pair of straight lines 
connecting the midpoint of the root chord to the tips. 
Additional computations of the sidewasli have been made 
using lifting lines composed of straight lines connecting the 
tips with the c/4 point and connecting the tips with the 3c/4 
point but the agreement with the exact sidewasli was not so 
good as that evidenced in figure 7. 

It is of interest that the span wise center of loading of the 
loading distributed along the lifting line connecting the c/2 
point with the tips (the lifting line yielding the best agree- 



Figure 7. — Comparison of lifting-line and lifting-surface sidewasli for 
0o=l.OO and 0.40. Circles represent points calculated by lifting- 
line method. 

ment with the lifting-surface results) was located longitudi- 
nally closer to the actual wing center of loading at 3c/4 than 
it was when the loading was distributed on the other two 
lifting lines. 

The wing loading in the examples just discussed was dis- 
tributed on one lifting line. Sidewasli obtained by distribut- 
ing the wing loading on more than one lifting line would 
probably show better agreement with the exact results in 
region E, because some effect of the longitudinal distribution 
of loading over the wing could then be realized. 

EXAMPLE OF FLOW-FIELD EFFECT ON VERTICAL TAIL 

The effect of the induced sidewasli velocity behind a 
rolling wing on the forces and moments contributed by a 
vertical tail can best be illustrated by analyzing a specific 
wing-tail configuration. The pertinent geometric char- 
acteristics of the wing-t ail model are (see fig. 8) : 
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Wing aspect ratio. _ 

Tail aspect ratio 



3. 2 



_ 1. 866 

Tail area 



__ 0.21 

Wing area 




Tail span 



0. 7 

Wing semispa 
Tail chord 



0. 6 

Wing chord 

Center-of-gravity location. _ 



_ . 0. 5c 


A free-stream Mach number of 1.6 (0=1.25) has been chosen. 
The wing leading edge for this Mach number is sonic (0o=l) 
and the vertical-tail leading edge is supersonic. Induced 
side-force and yawing moments for a number of longitudinal 
positions of the vertical tail have been determined by 
numerical integrations in a manner similar to that used in 
reference 11 to obtain the contribution of horizontal tails 
with supersonic leading edges to the lift and pitching moment. 
In making the numerical integrations, sidewash curves (figs. 
4 and 5) were used which had the infinity at the tip-cone 
intersection line faired through. Isolated vertical-tail forces 
and moments have been computed from the formulas given 
in reference 12. 

Figure 9 shows in stability-derivative form the variation 
of the induced, isolated, and total forces and moments with 
the longitudinal location of the vertical tail. For the 
example configuration chosen, the induced forces and 
moments are greater than the “isolated” forces and give rise 
to a positive C Yp coefficient and a negative C Up . Obviously, 
from the sidewash curves, if the vertical tail were moved 
away from the z-axis, the induced force would be reduced. 
The isolated forces and moments, on the other hand, w T ould 
increase and the total C Y and total C n would become 
negative and positive, respectively. 


ASSUMPTIONS AND LIMITATIONS 

In some cases, the assumptions made in the analysis, by 
necessity or for convenience, to allow the determination of 



Figure 8. — Wing-tail model used to illustrate effect of wing sidewash 
on a vertical tail. 



1.2 1.4 1.6 1.8 2.0 2.2 2.4 

Distance of tail apex behind wing apex measured in wing chords 


Figure 9. — Variation of induced, isolated, and total forces and 
moments (in stability-derivative form) with longitudinal location of 
vertical tail. 


the sidewash behind triangular wings by linearized super- 
sonic flow theory limit the application of the results. Some 
discussion of these assumptions and limitations may be 
useful. 

The validity of the assumption of a flat vortex sheet for 
wings with very low' aspect ratios is questionable, but, in 
the absence of experimental and theoretical information 
directly concerned with the vortex sheet behind rolling 
wfings, no definite statement can be made as to the effects 
that w T ing aspect ratio, roll velocity, and distance behind 
the trailing edge will have on the rolling-up of the vortex 
sheet. It may be possible, as suggested in reference 13, to 
get some indication of these effects from the data published 
in references 14 and 15 concerning the rolling-up of the vortex 
sheet behind wdngs at an angle of attack. 

The vortex sheet has been assumed not to rotate. The 
angle (in degrees) through which the vortex sheet w'ould 
rotate in moving from the wing trailing edge to a point d 
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distance downstream is given by 
. 180 pd 

180 pb / d \ 

7r 2V \b/2/ 

Substitution into this formula of values of j)b/2V usually 
encountered and of distance d up to two semispans will 
generally yield rotation angles small enough to be neglected. 

The sidewasli expressions derived in the analysis for rolling 
wings at zero angle of attack are also applicable for finite 
angles of attack because angle-of-attack loadings are sym- 
metrical and do not contribute to the sidewasli in the 22 - 
plane. The displacement of the assumed flat vortex sheet 
from its zero angle-of-attack position, however, must be 
accounted for; that is, the sidewasli given for a point x, z=0 + 
for the zero angle-of-attack case represents the sidewasli at 
the vortex sheet when the wing is at an angle of attack. At 
a distance d behind the trailing edge, the displacement of 
the vortex sheet below the trailing edge may be found 
(see sketch 3) from 



/ 

Wing chord extended 


Sketch 3 

Values of tan e for a triangular wing are given in reference 2. 


CONCLUDING REMARKS 

The variation of sidewasli with longitudinal distance in 
the vertical plane of symmetry behind rolling triangular 
wings traveling at supersonic speeds has been derived by 
linearized lifting-surface and lifting-line methods. The 
range of supersonic Mach numbers for which the lifting- 
surface results are valid is limited by the condition that the 
wing leading edges must be subsonic. The variations of 
lifting-surface sidewasli are presented in graphical form for 
a number of values of 0 O , a Mach number — leading-edge- 
sweep parameter. Sidewasli calculated from the lifting-line 
formula has been compared with the lifting-surface sidewasli 
for values of 0 O of 0.40 and 1.00. This comparison shows 
very good agreement of the lifting-line results with lifting- 
surface results except at the higher vertical distances close 
behind the trailing edge for 0 o =l.OO. The curves for 
0 o =O.4O reveal that, as 0 O is decreased from 1.00, the agree- 
ment close behind the trailing edge improves. 

An illustrative calculation of the sidewash-induced force 
on a half-delta vertical tail operating behind a triangular 
wing indicates that the induced force acts in opposition and 
is comparable in magnitude to the damping force created 
on the isolated rolling tail. In order to determine the total 
force which would act on the vertical tail of a steady-rolling 
aircraft in flight, an additional force induced by the aileron 
sidewasli should be calculated. This force may be of the 
same order of magnitude as, and opposed to, the force 
induced on the tail by the wing sidewasli. No attempt has 
been made in the present report to evaluate aileron sidewasli. 

Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., October 21 , 1955. 
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APPENDIX A 

TREATMENT OF SINGULARITIES IN v P , E AS GIVEN BY EQUATION (20) 

and a parts integration of expression (A2) gives 


In order to isolate the “infinite” terms, equation (23) may 
be written as 

© © 
zIIpV f' dxi , „ f r dx i E, 

x LJ.^ fe(-2K,+E,)+ J.iS*-T^ + 


V p, 


© © 

C c dx i JL ( _jr , FU f ‘ dx.\ E± ~| 

J/OoW 1 ,+ ,} J, fo* l-fc*J 


(Al) 


Terms ® and © are integrable by numerical methods, 
whereas terms © and © contain a first-order infinity at the 
limit x x =j (when x^f, k l =k 2 =l). 

Consider term © in equation (Al) 

f dx i k 2 E 2 


C f dx x k 2 l 

Jo 0o x i 1 


If k 2 is replaced by its equivalent (see symbols), term © 
becomes 


j: 


dx i 


(x—Xi) 2 —P 2 Z 2 —d 0 2 X l 2 


E 2 ^(x—Xi) 2 —[ Q 2 z 2 


(A 2) 


Assume expression (A2) to be of the form u dv , where 


j; 


and 


Then, 


u=Eo^{x-x l ) 2 -p 2 z 2 
dx i 


dv— 


v= 


-(x — Xi) 2 —P 2 Z 2 — d 0 2 X x 2 
1 ,1 — (1— floMi + S 


coth' 
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du= 


Mi 


E^jx-xtf-pz 2 ^ n1|> -(1— d 0 2 )xi+x 
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^do 2 x 2 +p 2 z 2 (l—do 2 ) r do(E 2 —K 2 ) K 2 k 2 (x— Xi) ~] A r 

Jo ^d 0 2 x 2 -\-fi 2 z 2 (l—d 0 2 ) L &2 Mi J 

(A3) 

Substituting for k x its equivalent in term © of equation 
(Al) results in 

6 0 XiEi dx x 


With 

and 


-// 


dv= 


doW—ix—x^+^z 2 
u=6 0 x 1 E l 
—dx i 


(A4) 


OqX 2 — (x—Xi) 2 -\-/3 2 z 2 
a parts integration of expression (A4) yields 

tanh - 


Mi E x _! - (1 — Qsi+s 


^9 o 2 x 2 +0 2 z 2 (l-6 o 2 )' " ^9 0 2 x 2 +p 2 z 2 (l-9 0 2 ) 
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tanh -1 


\e 0 2 x 2 +p 2 z 2 {i— e t 
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© L 1 0 0 O *A 2 J 1 
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The nonintegral term iii expression (A5) may be written in 
slightly different form as 


Eyy/^C—X i) 2 — j 8 2 3 2 tn|||| _, — ( 1 — 6£)Xi±X 

kx^e 0 2 x 2 +i3 2 z 2 a~do*) V<»o 2 x 2 +^2 2 (i-0o 2 ) 


(A6) 


Inasmuch as when Xi =/, k,=k 2 =l, it is clear that the evaluation of the integrated term in expression (A5) (equivalent 
to expression (A6)) at the limit x x —j will cancel the integrated term of expression (A3) evaluated at this same limit. 
The complete expression for r Pi is now seen to be 
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This equation in nondimensional form is presented in the analysis of the report as equation (24). 
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APPENDIX B 

EVALUATION OF vw.v 

The contribution of the doublets distributed in the wake to the sidewash in region D is given by equation (28) as 
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which by use of partial fractions can also be written 
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For ease in writing, — ^ has been replaced by a in equation (B2). Consider first the integral 


j:, 


\ 1 —s 2 ds 


(1 T"d 2 s 2 ) V 1- — ^* 2 , oV 

If the variable transformation s=sin 0 is made, expression (B3) becomes 

' T/2 cos 2 0 dd 


/; 


(1+d 2 sin 2 0) Vl— & 2 ,o 2 sin 2 0 
The evaluation of expression (B4) is given by formula (9), table 61 of reference 16 as 
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Expression (B4) could also have been integrated without recourse to tables with the aid of reference 17 (pp. 134 to 
136). The integration of expression (B3) 


Vl— s 2 ds 


(l+dV) 2 Vl-fe,o 2 s 2 

may be performed by using the relationship (see p. 13 of ref. 17 and p. 79 of ref. 18): 
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From equations (Bo) and (B6), equation (B7) becomes 
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Carrying out the differentiation in equation (B8) results in 
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(B7) 


(BS) 


G(k 2t0 ,d) + 

d V(d 2 + 1 ) (d 2 + k 2t o 2 ) J 


1 ^a 2 -\-l [~ ld 2 -\~k 2f o 2 K 2t() E 2 ,q 

2VaH^ 2 L 2 ’°V(a 2 +l) 3 V(d 2 +l)(d 2 +I 2 y)J 
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The remaining integration needed to evaluate v w ,d (cq. (B2)) is 

M Vl —s 2 ds 


X' 


(l-t 2 .oV ) 3/2 

Tli is integration may be reduced to standard elliptic forms by the transformation s=sn u and integrated to give 


x: 


1 tJ1—s 2 ds K 2t0 —E 2f0 
(i-k 2t0 vy' 2 ~ h, o 2 


(BIO) 


The sidewash v w ,d is now completely defined by equations (B2), (B5), (B9), and (BlO) as 


Vw, D~ 


3H(x—c)pk 


! t0 q f aVa 2 +l , 1 I/^/t n Vfl 2 1 [~ zr V a,2 ~l~^3, o 2 K 2t0 E 2t0 "j\ 

\L(^ 2 +^2 f o 2 ) 3/2 a ^ a 2 + 1 ^ a 2 + k 2t o 2 J 2 ’°’ V a2 +fe,o 2 L 2 *° ( a2 + 1 ) 3/2 Va 2 +1 V^ 2 +fe,o 2 J J 


2t 


H(x—c) z ak 2t o 3 4 5 6 7 8 9 f -Va 2 +1 , 3aV^ 2 +l 


a(a 2 +& 2f o 2 ) 3/2 2 (a 2 +^ 2 ,o 2 ) 5//2 2a^a 2 -\- 


1 (a 2 +fe.o 2 ) 3/2 ]' 


O(k 2t0 , a)- 


2(a 2 +I)(a : 


-K~2, 0 3 (g 2 , 0 o) *\ m,il 

L)(a 2 +^2,0 2 ) 2(a 2 +fe,o 2 ) 2 / u 


Considerable simplification of equation (Bll) may be accomplished by combining terms and noting that 


2 , , 2 ( x—c) 2 k 2 ,o 2 _ (x—c) 2 k 2 , 0 2 a 2 

*+*-2.0- p z 2 — do 2 c 2 


and 


V\v,d- 


II (3 6qc r 2-\ -a 


T [_a^a 2J rl 


G(k 2t o,a)- 


a(x—c)k 2t o 

Qqc(cl 2 -\-\) 


K, 


-] 


(B12) 


Replacing a by its nondimensional equivalent 


&Z Zq 


and H by its equivalent and then nondimensionalizing xby c and v w ,d by pb/ 2 gives 


V\v, 

pb / 2 


d 2 r22o 2- bl n ( 1 1 \ fcokUxo 1 )^n I 

/2 = ^w[vi7+l 2,0 ’W Wzo 2 + 1) J 


(B13) 


Equation (B13) is identical with equation (29) of the analysis, with the exception that, in the analysis, the 


function G 




has been written out. 
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